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Abstract 

i— i A semi-empirical three-dimensional model of turbulence in the approxima- 

tion of the far turbulent wake behind a body of revolution in a passive stratified 
i-^ medium is considered. The sought quantities are the kinetic turbulent energy, 

^ kinetic energy dissipation rate, averaged density defect and density fluctuation 

ypj variance. The full group of transformations admitted by this model is found. 

^ The model is reduced to the system of the ordinary differential equations due to 

O similarity presentations obtained and B-determining equations method. Sys- 

tem of ordinary differential equations satisfying natural boundary conditions 



C3 



was solved numerically. The solutions obtained agree with experimental data. 



1 Introduction 



CO 

> 

The turbulence play an important role in the formation of the ocean structure 
[TJ[2]. For example, the role of turbulence on the evolution of the spatial structure 
of a thin phytoplankton layer was examined in |3j. 

Semi-empirical models of turbulence are now widely used in methods of calcu- 
lation turbulent flows. However practically there are few analytical approaches to 
O research of this models. 

J> One of the examples of a three-dimensional free turbulent flow is a turbulent 

wake behind a body of revolution in a stratified medium. Sufficiently complete 
experimental data on the dynamics of a turbulent wake behind a body of revolution 
in a linearly stratified medium were obtained by Lin and Pao and presented in j2J. 

The turbulent wake behind an axisymmetrical body in a linearly stratified medium 
was numerically simulated in [5j. Based on hierarchy of semi-empirical turbulence 
models of second order, the numerical simulation of the dynamics of a turbulent wake 
in a stable stratified medium was carried out by Chernykh et al. |6]. A satisfactory 
agreement with experimental data (2) was obtained in (5j|6). 



A series of papers 11, 12, 16) was devoted to construction of similarity solutions 



of semi-empirical turbulence models. The present paper is a continuation of our 
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investigations. In this paper we consider three-dimensional semi-empirical model 
of the far turbulent wake behind an axisymmetric self-propelled body in a passive 
stratified medium |6~8||. Considering problem is equivalent to the problem of the 
development of a turbulent mixing zone in a passive stratified medium (7j[8]. 

In section 3 we have to define the admissible differential operators of the point 
groups of transformations 13,13] for considering model, which will allow us to pass 
to the system of the degenerate elliptic equations. In section 4 we focus on the 
solutions of the second order to the B-determining equation |15ij for degenerate 
elliptic equations. This gives the corresponding differential constraints and allow us 
to pass to the system of ordinary differential equations. In section 4 we will present 
the calculation results. 



2 Model 

To calculate the characteristics of the far turbulent wake behind an axisymmetric 
self-propelled body in a passive stratified medium we use the three-dimensional 
semi-empirical turbulence model |6j-[8] 
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In this equations uq is the velocity of an incoming undisturbed flow, e(x, y, z) is the 
turbulent kinetic energy, e(x, y, z) is the kinetic energy dissipation rate, (pi)(x, y, z) 
is the averaged density defect, and (p' 2 )(x,y,z) is the density fluctuation variance. 
The quantities C e = 0.136, C e = 0.105, C e2 = 1.92, C p = 0.208, C lp = 0.087, 
Ct = 1.25 are generally accepted empirical constants [9|[l0] . 

In what follows, we assume that the velocity of an incoming undisturbed flow 
equals unity. The marching variable x in equations Q-Q acts as the time. 

By analogy with 11,12 , for the model Q-Q we have to define the admissible 



differential operators of the point groups of transformations. 



3 Similarity solutions 

Group analysis of the system ([IJ-Q performed by a standard scheme 13,14 . 
The infinitesimal symmetry group of the model (JIJ-Q is spanned by eight vector 
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fields 
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Next consider the linear combination of the scaling vector fields Xq and Xj 
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The solution of the model ([TJj-Q invariant with respect to operator Z has the form 

e = x 2a - 2 E(£, rj), e = x 2a ~ 3 G(t, r,), ( Pl ) = x a H(£, r,), (p' 2 ) = x 2a R(£, rj), (5) 

where £ = y/x a , r\ = z/x a is the similarity variables. Substituting presentation ([5]) 
into Q-Q, we obtain the reduced system 
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Numerical analysis |6j-[8] of degeneration of the far turbulent wake in a passive 
stratified medium show that the functions E and G must be presented in the form 



E(t, V ) = E( v ^ 2 + V 2 ), G(t,v) = G( v ^ 2 + V 2 ) 



(10) 



Note that the presentations (10) are satisfied all the reduced equations (|6|)-([9|). 
Changing to polar coordinates £ = rcos(0), rj = r sin (</>), and by virtue of (10) the 
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reduced system become 
E 
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G' - 2E' sin(0) = 0, (13) 



AC P ^ (sm^Hr + C ^-H <t> ) = 0, (14) 



where E = E(r), G = G(r), H = H(r,<f>), R = R(r,<p). Here and elsewhere, 
subscripts denote derivatives, so H r = dH/dr, etc. We now apply the BDE method 



15 to reduce the equations (13), (14) to some ordinary differential equations. 



4 BDE method 



Consider more general equation than (13) 

H M + r 2 H rr + A(r)H r + B(r)H + C{r) sin(0) = 0, 



(15) 



where A(r), B(r), C{r) are arbitrary functions. We take the B-determining equation 



corresponding to (15) of the form 

D 2 h + r 2 D 2 h + 6i(r, <j>)D r h + 6 2 (r, <j>)h = 0. 



(16) 



Here and throughout D^, D r are the operators of total differentiation with respect 
to <p and r. The functions b\(r, <fi) and b 2 (r, <ft) are to be determined together with 
the function h. Note that for classical determining equations 13,14] holds 



&i(r,0) = A(r), & 2 (r,. 



B(r). 



We seek second order solution of (16) of the form 

h = Haa, + hi U, H, Ha 



(17) 



Substituting (17) into BDE (16) leads to an equation which includes derivatives of 



the fourth order. We can express the derivatives H rr AA, H^a, H^, by 



means of (15). Setting the coefficient of H rrr equal to zero we obtain bi(r, 0) = A(r). 
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The left-hand side of (16) is a polynomial with respect to H rr and H r &. This 



polynomial must identically vanish. Collecting similar terms we obtain the equations 

= 0, h lHH ^ = 0, 2 (A(r)H r + B(r)H + C(r) sin(0)) h lllAH - 

2H <j) h lHH 



2h 1 ^+B(r)-b 2 (r,cj ) ) = 0. 



It is easy to show that the general solution of the equations (18) is 

h (<p, H, H+) = h^H* + h 3 (<f>, H), b 2 (r, <f)) = B(r) - 2h' 2 {(j>). 



Substituting the functions 6 1( b 2 an d h\ into BDE (16) we obtain that the left 



hand side of (16) is a polynomial with respect to H r and H^. This polynomial must 



identically vanish. Collecting similar terms leads to the following equations 

h 3HH = 0, 2h^ H + K((f>) - 2ti 2 {<P)h 2 {<j>) = 0, (B(r)H + C(r) sm(<l>))h 3B - 

+ (2h! 2 (cf>) - B(r))h 3 + C(r)(cos(0)/i 2 (0) - sin(0)) = 0. 

The equations (19) imply 

h 3 (<j>,H) = Q^ 2 (0) 2 - hi'^) + H, 
ti 2 {<j>) - h 2 2 (<p) - 2 cot(0)/i 2 (0) + 2(1 - h A ) = 0, 
here /14 is arbitrary constant. 



(19) 



Clearly, that the Riccati equation (20) has the partial solution 



h 2 (4>) = tan(0) 

for h 4 = 1/2. 

Thus we find the second order solution of the BDE (16) 

h = H H + tan(0)^. 
The corresponding differential constraint h = has the general solution 

H = H 1 (r)sin(4>)+H 2 (r), 
where Hi and H 2 are arbitrary functions. 
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(21) 
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Substitution (22) into equation (14) gives 
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By analogy with the case of the equation (13), consider more general equation 



than (23) 



R M + r 2 R rr + K(r)R r + L(r)R + M(r) sin 2 (0) + N(r) sin(0) + P(r) = 0, (24) 



where K(r), L(r), M(r), N(r) and P(r) are arbitrary functions. The BDE method 



applied to equation (24) gives rise to the following results: 
bi (r, 4>) = K(r), b 2 (r, <p) = L(r) - 



sin 2 (2(f)) 



h = - 2cot(24>)R 4> , 
N(r) = 0. 



(25) 
(26) 



The formula (26) for the equation (23) takes the form 

{H[ - l)H' 2 = 0. 
Clearly, that we must explain the case 

H' 2 = 0. 



(27) 



Integrating differential constraint h = corresponding to the BDE solution (25), we 
find 

R = Ri(r) sin 2 (<f>) + R 2 (r), (28) 
where R\(r) and R 2 (r) are arbitrary functions. 



Thus in the similarity variables £ and r\ from (22), (|27[), (28) we have 



H(^ v ) = H 3 (y/e + V 2 )v + H 2 , R& rj) = R 3 (Ve + V 2 W + R2(Ve + V 2 ), (29) 



where H 3 = H x /^i 2 + r} 2 and R 3 = R 1 /(£ 2 + r/ 2 ). 

This allow us to reduce the model (jlj)-( |4"|) t o the system of ordinary differential 
equations. Substituting presentations (10), (|29|) into the reduced system ([6])-([9]) we 
obtain 



H 2 = 


0. 


E" = 


E> 


G" = 


a 


m= 


H' 3 




R's 


2C P H' 3 ( 






R 2 = 


R 2 



G' E' 
~G~ 2 E~ 



E 



G 

G E 
G' 



1 

T 

E' 1 

T 

E' 3 

T 

E' 5 



G 
~C^E 
G 
' C t E 2 
arG 

~C P ~E 2 
arG 



G E' 
2(a - 1) + - - ar— 



(2a-3)G + 



C e2 G 2 
E 



— arG' 



G E 
[Hs ~ 1) 

T 



r C l0 E 2 



H 3 - 




T 




2R 3 


/ G' 


T 


\G 



G 



E' 



(30) 
(31) 
(32) 



E' C t tG 2 
2^ + 



E 2C lp E 3 



m 



3 ' 



(33) 



arG 



a e' i 

~G ~ 2 ¥ ~ t ~ C l0 E 2 



+ 



RoG ( CtG 



C l0 E 2 V E 



+ 2a - 



2R 3 



2C n 



where r = a/£ 2 + rf 



(Hs - If 



(34) 



6 



5 Calculation results 



System (30)-(34) has to satisfy the conditions 

E' = G' = H[= R[ = R' 2 = 0,t = 0, 
E = G = E x = Rx = R 2 = 0, r oo. 



(35) 
(36) 



Conditions (35) takes into account flow symmetry with respect to the OX axis. 



The boundary conditions (36) imply that all functions take zero values outside the 
turbulent wake. 





(a) 



(b) 





(c) (d) 

Figure 1: Calculated profiles as £ = 0: (a) normed profile of E, (b) normed profile 
of G, (c) profile of H, (d) normed profile of R. 



The system (30)-(34) of ordinary differential equations satisfying boundary con- 



dition (35), (36) was solved numerically. Additional difficulties are caused by the 
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fact that the coefficients of ordinary differential equations have singularities. The 
problem was solved by a modified shooting method and asymptotical expansion of 
the solution in the vicinity of the singular point |16|. 




Value of a a taken to be equal to 0.23. The results for the problem solu- 
tion are illustrated in Figs. 1, 2. Figure 1 shows the profiles of the functions 
E/E ,G/Gq, H and R/Rq as £ = 0, where subscript denote axial value. The 
functions E/E ,G/G , H and R/Rq are plotted in Fig. 2. 

The function H(0,i]) characterizing the degree of fluid mixing in the turbulent 
wake a given in Fig. lc. As can be seen, the maximum value of this function slightly 
differ from 0.25, which is consistent with the present notions of incomplete fluid 
mixing in the wakes \L7\. 
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Figure 3: Axial values of the turbulent energy. 

In Fig. 3 the axial values of the turbulent energy are compared with Lin and Pao's 
experimental data |4j, Hassid' s computational results |5] and results of numerical 
computations (6j[7]. We have borrowed this figure from work [6] and have put the 
values. We can see satisfactory agreement with Lin and Pao's experimental data 
here as well. 

Conclusion 

The main results of the paper are as follows. The three-dimensional semi- 
empirical turbulence model of the far turbulent wake behind an axisymmetric self- 
propelled body in a passive stratified medium was reduced to the system of ordinary 
differential equations due to similarity presentations obtained and B-determining 
equations method. The system of ordinary differential equations satisfying natural 
boundary conditions was solved numerically. The solutions constructed agree with 
experimental data. 
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